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Abstract 



We extract the ratio of the down (d) and up (u) parton distribution functions (PDF's) 
from the ratio of NMC deuteron and proton structure function F2 /F2, using corrections 
for nuclear binding effects in the deuteron, which are extracted from the nuclear depen- 
dence of SLAC F2 data. Significant corrections to the d quark distribution in standard 
PDF's are required, especially at high x. The corrected d/u ratio is in agreement with the 
QCD prediction of 0.2 at x = 1. The predictions for the recent CDF W asymmetry data 
using PDF's with the corrected d/u ratio give much better agreement at large rapidity. 
Using the updated d/u ratio and the most recent world average for a^, we perform a NLO 
global fit to all DIS data for F2 and R, and estimate the size of the higher twist contri- 
butions using both a renormalon model and an empirical model. We find that with the 
updated value of as , the magnitude of the higher twist terms is half the value of previous 
analysis. With the inclusion of target mass and higher twist corrections, the standard 
NLO PDF's with the updated d/u ratio describe the SLAC F2 data up to x = 1.0. When 
the analysis is repeated in NNLO, we find that the additional NNLO contributions to R 
account for most of the higher twist effects extracted in the NLO fit. The analysis in 
NNLO indicates that the higher twist effects in R, F2 and xF^ (e.g. GLS sum rule) are 
very small. 



* To be published in proceedings of the 6th International Workshop on Deep Inelastic Scattering and QCD, 
Brussels, Apr. 1998. Email: ukyang@fnal.gov 
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Figure 1: [a] The total correction for nuclear effects (binding and Fermi motion) in the deuteron, F^jF^^^ ■, 
as a function of x, extracted from fits to the nuclear dependence of SLAC F^ electron scattering data, 
[b] Comparison of NMC Fl^'^'^ jF'^ (corrected for nuclear effects) and the prediction in NLO using the MRS(R2) 
PDF with and without our proposed modification to the dju ratio. 



1 Introduction 

Recent work on parton distributions functions (PDF's) in the nucleon has focussed on probing 
the sea and gluon distribution at small x. The valence quark distributions has been thought to 
be relatively well understood. However, the precise knowledge of the u and d quark distribution 
at high X is very important at collider energies in searches for signals for new physics at high . 
In addition, the value of (i/w as x ^ 1 is of theoretical interest. Recently, a proposed CTEQ 
toy model El included the possibility of an additional contribution to the u quark distribution 
(beyond x > 0.75) as an explanation for both the initial HERA high anomalyi, and for the 
jet excess at high-Pj at CDPi. In this communication we conclude that a re-analysis of data 
from NMC and SLAC leads to a great improvement in our knowledge of the PDF's at large x. 



2 Extraction of d/u at high x 

Information about valence quarks originates from proton and neutron structure function data. 
The u valence quark distribution at high x is relatively well constrained by the proton structure 
function • However, the neutron structure function F^, which is sensitive to the d valence 
quark at high x, is actually extracted from deuteron data. Therefore, there is an uncertainty 
in the d valence quark distribution from the corrections for nuclear binding effects in the 
deuteron. In past extractions of from deuteron data, only Permi motion corrections were 
considered, and other binding effects were assumed to be negligible. Recently, the corrections 
for nuclear binding effects in the deuteron, F2 / F2^^ ^ have been extracted empirically from fits 
to the nuclear dependence of electron scattering data from SLAC experiments E139/140i. The 
empirical extraction uses a model proposed by Prankfurt and Strikmani, in which all binding 
effects in the deuteron and heavy nuclear targets are assumed to scale with the nuclear density. 
The correction extracted in this empirical way is also in agreement (for x < 0.75) with recent 
purely theoretical calculations! of nuclear binding effects in the deuteron. The total correction 
for nuclear binding in the deuteron is about 4% at x = 0.7 (shown in Pig. |l|(a)), and in a 
direction which is opposite to what is expected from the previous models which only included 
the Permi motion effects. 

The ratio F^/F^ is directly related to d/u. In leading order QCD, 2F^/F^ - 1 ~ 
(1 + Ad/u) /(4 + d/u) at high x. We perform a NLO analysis on the precise NMC Fi/F^ data§ 
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Figure 2: [a] The d/u distribution at Q^=1Q as a function of x. The standard MRS(R2) is compared to our 
modified MRS(R2). [b] Comparison of the CDF W asymmetry data with NLO standard CTEQ3M, MRS(R2), 
and modified MRS(R2) as a function of the lepton rapidity. The standard CTEQ3M with a resummation 
calculation is also shown for comparison. 



to extract d/w as a function of x. We extract the ratio "/ -Ff by applying the nuclear binding 
correction F^/F^^'p to the F^/F^ data. 

As shown in Fig. |l](b), the standard PDF's do not describe the extracted Ff^'^/Fl'. 
Since the u distribution is relatively well constrained, we find a correction term to d/u in the 
standard PDF's (as a function of x) by only varying the d distribution to fit the data. The 
correction term is parametrized as a simple quadratic form, 6{d/u) = (0.1 ± 0.01) (x + l)a; 
for the MRS(R2) PDF, where the corrected d/u ratio is {d/u)' = {d/u) + 6{d/u). Based on 
this correction, we obtain a MRS(R2)-modified PDF as shown in Fig^(a). The corrections to 
other PDF's such as CTEQ3M is similar. The NMC data, when corrected for nuclear binding 
effects in the deuteron, clearly indicate that d/u in the standard PDF's ii is significantly 
underestimated at high x as shown in Fig. 0. Fig.|^(a) also shows that the modified d/u ratio 
approaches 0.2 ± 0.02 as x — 1, in agreement with a QCD prediction 0. In contract, if the 
deuteron data is only corrected for Fermi motion effects (as was mistakenly done in the past) 
both the d/u from data and the d/u in the standard PDF's fits approach as x ^ 1. 

Information (which is not affected by the corrections for nuclear effects in the deuteron) 
on d/u can be extracted from W production data in hadron colliders. Fig.H(b) shows that the 
predicted W asymmetry calculated with the DYRAD NLO QCD program using our modified 
PDF is in much better agreement with recent CDF data0 at large rapidity than standard 
PDF's. 

When the modified PDF at Q^=1Q is evolved to Q^=10, 000 using the DGLAP NLO 
equations, we find that the modified d distribution at x = 0.5 is increased by about 40 % in 
comparison to the standard d distribution. The modified PDF's have a significant impact on 
the prediction of the cross sections in the HERA high region because the charged current 
scattering with positrons is on d quarks only. It also impacts neutral current scattering because 
of the large coupling of the Z to d quarks at high Q^. The modified PDF's also lead to an 
increase of 10% in the QCD prediction for the production rate of very high jets in hadron 
colliders. 

3 Higher twists effects at high x 

Since all the standard PDF's, including our modified version, are fit to data with x less than 
0.75, we now investigate the validity of the modified MRS(R2) at very high x by comparing 
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Table 1: Results of the higher twists fits (with pQCD+TM) to the DIS F2 data 





SLAC(%) BCDMS(%) a b c 


Proton 
Deuteron 
BCDMS main sys (A) 


-1.2 -4.1 0.50 3.2 0.11 
-0.6 -2.0 0.35 1.5 0.26 
1.35 



to F2 data at SLAC. Although the SLAG data at very high x are at reasonable values of 
(7 < < 31 GeV"^), there are in a region in which non-perturbative effects such as target 
mass and higher twist are very large. We use the Georgi-Politzer calculation for the target 
mass corrections (TM). These involve using the scahng variable ^ = 2x/(l + \Jl + AM'^x'^ / Q"^) 
instead of x. Since a complete calculation of higher twist effects is not available, the very low 
data is used to obtain information on the size of these terms. 

We use two approaches in our investigation of the higher twist contributions: an em- 
pirical method, and the renormalon model. In the empirical approach, the higher twist con- 
tribution is evaluated by adding a term h{x)/Q'^ to the perturbative QCD (pQCD) prediction 
of the structure function (including target mass effects). The x dependence of the higher 
twist coefficients h{x) is fitted to the global DIS F2 (SLAC, BCDMS, and NMC) datafflS 
in the kinematic region (0.1 < a; < 0.75, 1.25 < < 260 GeV"^) with the following form, 
F2 = ^P'3'^-^+'^^(x _|_ h{x) / Q"^) f (x) . Here f{x) is a fioating factor to investigate possible x de- 
pendent corrections to our modified PDF. A functional form, 0(^3^ — c) for h{x) is used in the 
higher twist fit to estimate the size of the higher twist terms above x = 0.75. The SLAC and 
BCDMS data are normalized to the NMC data. In the case of the BCDMS data, a systematic 
error shift A (in standard deviation units) is allowed to account for the correlated point-to-point 
systematic errors. The empirical higher twist fits with the modified NLO MRS(R2) pQCD pre- 
diction with TM have been performed simultaneously on the proton and deuteron F2 data with 
11 free parameters (2 relative normalizations and 3 parameters for C{x) per target and the 
BCDMS A) . The results of the fit are given in table |I[ 

We find that empirical higher twist fit describes the data well {x^/DOF = 843/805). 
The size of the higher twist contributions in the proton and deuteron are similar. The magnitude 
is almost half of those extracted in previous analysis of SLAC/BCDMS dataS This is because 
that analysis was based on q;^(M|) = 0.113, while the MRS(R2) PDF uses as(M|) = 0.120, 
which is close to the current world average. In the renormalon model approach§: the model 
predicts the complete x dependence of the higher twist contributions to F2, 2xFi, and xF^, 
with only two unknown parameters A2 and ^4. We extract the A2 and parameters, which 
determine the overall level of the l/Q"^ and 1/Q^ terms by fitting to the global data set for F2 
and R{= ^2(1 -|-4Ma;^/ Q^)/2xFi — 1). The values of A2 and for the proton and deuteron are 
same in this model. The x dependence of 2xFi differs from that of F2 but is same as that of xF^ 
within a power correction of 1/Q^. Our fits can also be used to estimate the size of the higher 
twist effects in xF^ (e.g. the GLS sum rule). The higher twist fit in this approach has employed 
the same procedure as the empirical method. Fig. ^ shows that the model yields description 
of X dependence of higher twist terms in both F2 and R with just the two free parameters 
ix'/DOF = 1577/1045). The CCFR neutrino data0 is shown for comparison though it is 
not used in the fit. The extracted values of A2 are —0.093 ± 0.005 and —0.101 ± 0.005, for 
proton and deuteron, respectively. The contribution of A4 is found to be negligible. We find 
that the fioating factor f{x) for the deuteron deviates from 1 and is also bigger than that for 
the proton, unless the modified MRS(R2) PDF is used. This refiects our earlier conclusion that 
the standard d distribution is underestimated at high x region. As expected, the extracted 
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Figure 3: The description of higher twist fit using the renormalon model with the modified NLO MRS(R2) PDF. 
The CCFR neutrino data is also shown for comparison, [a] Comparison of F2 and NLO prediction with and 
without higher twist contributions, [b] Comparison of R and NLO prediction with and without the renormalon 
higher twist contributions. 



A2 value is half of the previous estimated value El of A2 based on SLAC/BCDMS (as(M|) 
= 0.113) analysis. Since both of these approaches yield a reasonable description for the higher 
twist effects, we proceed to compare the predictions of the modified PDF's (including target 
mass and renormalon higher twist corrections) to the SLAC proton F2 data at very high x 
(0.7 < a; < 1). 



4 Parton distributions functions at very high x 

There is a wealth of SLAC datai in the region up to x = 0.98 and intermediate {7 < < 
31 GeV'^). Previous PDF fits have not used these data. We use the estimate of the higher 
twist effects from the models, based on the data (below x < 0.75) described above. Note that 
the data for x > 0.75 is in the DIS region, and the data for x > 0.9 is the resonance region. It 
is worthwhile to investigate the resonance region also because from duality arguments o it is 
expected that the average behaviour of the resonances and elastic peak should follows the DIS 
scaling limit curve. Fig. ^ shows the ratio of the SLAC data at very high x to the predictions 
of the modified MRS(R2). With the inclusion of target mass and the renormalon higher twist 
effects, the very high x data from SLAC is remarkably well described by the modified MRS(R2) 
up to X = 0.98. The good description of the data by the modified MRS(R2) is also achieved 
using the empirical estimate {h{x) / Q"^) of higher twist effects as shown (dashed line) in Fig. ^(c). 
The data at x > 0.9 and 7 < < 11 GeV^ is somewhat lower than the predictions of 
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the modified MRS(R2) but relatively the large data {21 < < 30 GeV^) show a good 
agreement as shown in Fig. ^(b). This could be due to the missing elastic contribution as 
expected from duality arguments. Fig. ^(b) also shows that the CTEQ Toy model (with an 
additional 0.5% component of u quarks beyond x > 0.75) overestimates the SLAC data by a 
factor of three at x = 0.9 (DIS region). From these comparison, we find that the SLAC F2 data 
do not support the CTEQ Toy model which proposed an additional u quark contribution at 
high explanation of the initial HERA high anomaly and the CDF high-P^ jet excess. 

5 Conclusion 

We find that nuclear binding effects in the deuteron play a significant role in our understanding 
of d/u at high x. The modified PDF's with our d/u correction are in good agreement with the 
prediction of QCD at x = 1, and with the CDF W asymmetry data. With the inclusion of 
target mass and higher twist corrections, the modified PDF's describe all DIS data not only up 
to the very highest x, but also all the way down to Q"^ = 1 GeV^. 

Note that when the analysis is repeated in NNLOill, We find that the additional NNLO 
contributions to R account for most of the higher twist effects extracted in the NLO fit. The 
analysis in NNLO indicates that the highest twist corrections to the GLS sum rule are very 
small (the fractional contribution to the pQCD GLS sum rule is -0.0058/(5^ _ 0.013/Q^). 
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Figure 4: The ratio of SLAC data at high x to the predictions of the modified MRS(R2). The data 
are compared to the NLO pQCD prediction [a] with the additional target mass effects [b] and the inclusion 
of target mass and renormalon higher twist effects [c]. The dashed line is the ratio of the prediction with 
empirical higher twist effects to that with the renormalon higher twist effects. 
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Figure 5: Comparison of SLAC Ff data with the predictions of the modified MRS(R2), CTEQ4M and the 
CTEQ toy model at high x and higher (20 < < 3f). 
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